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SUMMARY 


Flight measurements of the dynamic longitudinal— stability cba.ra.cter— 
istics of a 35° swept-wing airplane are presented in this report. Pitch- 
ing velocity responses to pulse-type elevator motions were recorded over 
a Mach number range of 0.60 to 1.04 at an altitude of approximately 
35,000 feet. The measured period and damping of these oscillatory 

responses were used to compute the static stability parameter C„ , the 

% 

factor Cj^ + Cj,^, and the number of cycles for the oscillations to damp 

to l/lO amplitude. These results are compared with the same character- 
istics estimated from static wind-tunnel data. 

The test results show a sharp reduction in damping at a Mach number 
of 0.92 with maximum damping occurring at Mach numbers of 0.88 and O.^k. 
Above 0.9^ there is again a decrease in damping that continues over the 
remainder of the test range. The static stability increases rapidly with 
Mach number through the transonic range, except for a slight decrease at 
M = 0.92. 


INTRODUCTION 


At the present time there are only limited data available on the 
dynamic behavior of swept-wing aircraft, especially at transonic speeds. 
The flight-test program being conducted by the NACA on a 35° swept-wing 
airplane includes the determination of dynamic directional—, lateral—, 
and longitudinal— stability characteristics as affected by Mach number. 
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This report presents the preliminary results of a flight investiga- 
tion into the longitudinal-stability characteristics of the test airplane. 
The simple method of analysis used here is possible because the airplane 
exhibits a lightly damped short— period oscillatory response. The com- 
puted flight-test characteristics are compared with estimations made from 
static wind-tunnel data. 

More complete analyses, using additional flight data, are being car- 
ried out in order to determine the frequency responses of the airplane 
and the individual stability derivatives through the subsonic and tran- 
sonic speed ranges. 


SYMBOLS 


Cp lift coefficient 

C m pitching-moment coefficient about airplane center of gravity 

number of cycles required for oscillation to damp to l/lO ampli- 



tude 


I 


y 


moment of inertia about Y axis, slug feet squared 


M 


Mach number 


P 


period of oscillation, seconds 


S 


wing area, square feet 


S.|. horizontal— tail area, square feet 

V velocity, feet per second 


damping coefficient 


b w wing span, feet 



c 


wing mean aerodynamic chord 


, feet 


c 


local wing chord, feet 


i angle of incidence of horizontal tail, radians 
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k aerodynamic restoring moment 

Z t tail length, feet 


m 

q 

t 

y 

a 
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mass, slugs 

pitching velocity, radians per second 
time, seconds 

distance from plane of symmetry, normal to plane of symmetry, feet 

angle of attack, radians 

elevator deflection, radians 

downwash angle, radians 

tail efficiency 

atmospheric density, slugs per cubic foot 

dynamic pressure, pounds per square foot 

da 

dt 


C dC L 


U T 

L a 

da 

c 

dC^ 

Sx 

da 

c 

dc m 
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Subscript 


horizontal tail 


TEST EQUIPMENT 


The test airplane was a standard North American F— 86 a— 5 with exter- 
nal booms added as shown in figures 1 and 2. The pertinent physical 
characteristics of the airplane are listed in table I. 

Standard NACA optical instruments synchronized at l/lO— second inter- 
vals by a common timer were used to determine static and impact pressures 
on a Kollsman airspeed head mounted in the nose boom (fig. 2), normal 
acceleration at the center of gravity, pitching velocity, and elevator 
and stabilizer positions. The true Mach number and pressure altitude 
were obtained by correcting the static and impact pressures in accordance 
with the calibration of reference 1. 


TEST PROCEDURE AND REDUCTION OF DATA 


The flight program consisted entirely of elevator pulse-type maneu- 
vers and covered a Mach number range of 0.60 to 1.04 with corresponding 
variations in trim lift coefficient of 0.35 to 0.12. The airplane was 
trimmed in either level or diving flight at the start of each run and, 
after application of the pulse, the control stick was returned approxi- 
mately to the trim position and held against a chain stop so as to mini- 
mize any inadvertent elevator motion. The success of this system is 
shown in the typical flight records presented in figure 3. The pitching- 
velocity responses were analyzed by the method shown in appendix A. 
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All records at speeds greater than 0*92 Mach number were taken in 
diving flight. For the calculations presented herein, however, the alti- 
tude and Mach number were assumed to remain constant at their average 
values during the 2- to 4-^econd duration of each run. The validity of 
this assumption is borne out by the data shown in figure 3(a) wherein 
the Mach number does not vary over 0.7 percent during the interval in 
which the data were taken. The change in flight— path angle associated 
with a dive may have some effect on the response of the airplane; how- 
ever, analytic investigation showed that in this case the effect on the 
period and damping of the oscillatory mode was negligible. 

Since the actual flights covered an altitude range of 30,000 to 
^0, 000 feet, the period and damping measurements were corrected to a 
standard altitude of 35*000 feet by assuming that the damping varies 
directly with the air density and that the period is inversely propor- 
tional to the square root of the density. 


RESULTS AND DISCUSSION 


The period of the oscillations as measured from flight records is 
plotted in figure 4 as a function of Mach number. It should be remem- 
bered that all the data presented herein were obtained in either level 
or diving flight, the lift coefficient varying between 0.35 and 0.12 at 
35*000 feet over the speed range covered. It can be noted that there is 
a continuous decrease in period with increasing speed except for the 
slight rise between Mach numbers of 0.92 and 0.94. Also shown are values 
of period estimated from wind-tunnel data of reference 2 by the method 
shown in appendix B. The values computed from the wind-tunnel data are 
somewhat smaller in magnitude but show the same variation with speed as 
the test values of period. 

Figure 5 shows a plot of damping coefficients that were obtained 
from the oscillation records by the method described in appendix A. 
Examination of this curve shows that, after reaching a peak value at a 
Mach number of 0.88, the damping drops sharply to a low value at 0.92, 
after which there is a sudden increase to a maximum at a Mach number of 
0.94. Above 0.94 the damping decreases in an irregular fashion over the 
remainder of the test range. It should be noted that the sharp decrease 
•in damping is defined by five separate test points taken independently 
at about Mach number 0.92. This is also the same speed at which the 
irregularity in the period was noted. 

Also plotted in the same figure are estimated values of damping 
(see appendix B) which show somewhat smaller magnitudes over the range 
of Mach numbers covered by the wind tunnel. This difference nay result 
from the assumptions made regarding both de/da (0.5) and the portion 
of Mq due to the wing and fuselage (0.25 Mq^). 
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Considerable scatter can be noted in the test points of figure 5* 

This is an indication of the difficulty involved in fairing an exponen- 
tial envelope to the test data. The accuracy of the measurements is 
poorer at the lower speeds because the oscillations were almost completely 
damped within one cycle. It was very difficult to determine the damping 
from these test measurements at Mach numbers below 0.70. At the higher 
speeds (above Mach number 1.0) the scatter is not due entirely to inaccu- 
rate measurements since the records consisted of two or more well-defined 
cycles. The period measurements appear to be more consistent and no dif- 
ficulty was encountered except at the lowest test speeds. 

The period and damping test data (figs. 4 and 5) were used to comr- 
pute the static stability parameter shown in figure' 6. The flight— 

test curve shows a rapid increase in static stability with increasing 
Mach number except for a dip, which is similar to that previously noted 
in figure 4. 


The computed values of C m should have the same degree of accu— 

racy as the original period measurements since the period of oscillation 
depends almost entirely on C^. As indicated in appendix A, however, 

a small error was introduced by the omission of the term Z a M q /mVT y ; 
but in this case the error was estimated to be negligible at high speeds 
and less than 4 percent at the lowest test speed. 


Also plotted in figure 6 are values of obtained from the data 

of reference 2. There is fair agreement between the two curves, although 
the wind-tunnel results do show higher values of static stability corre- 
sponding to a shift in neutral point of about 2.5 percent. 


The variation of with Mach number is plotted in figure 7. 

This parameter was computed by combining experimental values of damping 
coefficient b taken from the faired curve of figure 5 with values of 
taken from figure 8. 1 These calculati ons are explained in appen- 
dix A. 


By comparing figures 5 and 7, it can be seen that both curves exhibit 
the same general characteristics at Mach numbers greater than 0.92. At 
the lower speeds, however, the influence of causes the total damp- 

ing b to increase with Mach number, while the term C m ^ + Cj^ is actu- 
ally decreasing. The sudden drop in b at Mach number 0.88 is due also, 

1 Figure 8 is the result of estimations made by the North American Avia- 
tion Company based on wind-tunnel (reference 2) and transonic— bump 
tests . 
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at least partially, to the decrease in Cr that occurs at this same 
speed. 


estimated 


Also plotted in figure 7 is the variation of 
from reference 2. Since the same values of C T were used in computing 

vX 

both curves in this figure , the poor agreement between the two is the 
result of discrepancies between the estinated and the experimental damp- 
ing coefficients shown in figure 5* In general, the accuracy to which 
C mq + C m^ can determined experimentally is limited by the accuracies 

of both b and C-^. As shown in appendix A, C^ + is propor- 

tional to the difference between b and Z a /mV and thus will usually 
be less accurate than either of these two quantities. Examination of the 
wind-tunnel data used in this report shows that Z^/mV is approximately 
one— half as large as b; thus a 1— percent error in b results in a 
2— percent error in + C^. 


In figure 9 is plotted the number of cycles required for the oscil- 
lations to damp to l/lO amplitude. The curve shows that below a Mach 
number of O.85, approximately O.95 cycle is required. Above this speed 
the number of cycles required increases rapidly with Mach number, exhib- 
iting irregularities similar to those noted in the damping curve. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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APPENDIX A 


METHOD OF AMLYSIS 


For the analysis used in this report, the following longitudinal 
equations of motion are assumed: 


mV (a — q) = Z a a + Z55 


(Al) 


I y q = MqCC + M&a + M q q + M§5 

Solving the two equations simultaneously for q/S yields 

q CjD + Cq 

5 D 2 + hD + k 


(A2) 


(A3) 


where 


c o “ 


Ci = 


dt 


(a 4) 

ZsMa " 

M 5 z a 

(A5) 

mV ] 

[ y 


Z&Md + 

% 

(A6) 

mVTy 

I y 


_Za_ 

Mn Ma 

—2 — = damping coefficient 

(AT) 

mV 

x y % 


z a M q _ 

M 

— = aerodynamic restoring moment 

(A8) 

mVIy 

I y 
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The damping coefficient b can be obtained from a typical oscilla- 
tion record such as shown in figure 10 where the equation of the envelope 
is assumed to be of the form 


a 



(A9) 


Thus 


b = ln(a-i/a g ) 
2 ~tz > — ‘*fci 


(A10) 


The aerodynamic restoring moment k is determined from the period 
the oscillation as follows: 


P 


2k 



and, by rearranging terms. 



(All) 


(A12) 


The term Cjj^ can be determined from k by using equation (A8), 

omitting the term Z a Mq/mVIy. Investigation has shown that, in this case, 
this term is very small compared to Ma/l . Thus, 


k 


% _ Sc pV 2 

I y I y 2 


(A13) 


and 

c - kI y 2 
- - sT 


Uik) 


The damping factor + Cj^ can be obtained from equation (A7) 

using known values of b and C T . 

h <x 
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-i- (M n + Mx ) = - — - b 
I ' S u mY 


(A15) 


Cm q + C % pVSc 2 2m 


UX 


7 ( Cj. 


(A16) 


The term Ci/ 1Q is obtained from the relation 


r. . _ In 10 

Ol/lO V- 

p — 

2 


(A17) 
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APPENDIX B 

ESTIMATION OF DYMMIC STABILITY PARAMETEBS 
FROM STATIC WIND-TUNNEL TESTS 


In. general, it is not possible to obtain dynamic stability parameters 
directly from static wind-tunnel tests. By making certain assumptions, 
however, it was possible to estimate the values of the coefficients b 
and k using the data of reference 2; then the other stability parameters 
were computed by the methods of appendix A. 

In appendix A it was shown that 


* Za M % 

k = 

mY! y I y 


To determine b and k it was necessary to evaluate 7 ^, Mq, 

M&, and Ma usin S wind-tunnel plots of Cj^, dC m /di t , and dC m /dC L . 
This was done as follows: 


For Z a , 


Z a 


C I<x 


pf_ 

2 


S 


For Mq, the portion of Mq due to the tail alone can he expressed as 




dCLt, PV 
da t 2 


dC m 


dC 


and since 


’^t ^t^t 
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then 

dC m pVScZ-^ 

M< Jt = iq — — 

In addition, however, the portion of Mq due to the wing and fuselage 
must also he estimated; in this case, it was assumed to he 25 — percent of 
Mq^.. Therefore, 


M 


q 


1.25 


dC m pVSc Z-t 
lit 2 


For M^, can he expressed as 


a ^ dot l a 2 


or 


M * = M n 
a q t da 


The term de/da was estimated to he 0.5 over the wind-tunnel speed range. 


For can he obtained from the plot of dC m /dC L in reference 2 

where C m is the moment coefficient about the quarter— chord point. The 

center of gravity of the test airplane was located at approximately 22.5 
percent of the mean aerodynamic chord. Thus, 


C 


“a 



VdC L 


0.025 


and 


Hx = c ma |^ Sc 
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TABLE I.- PHYSICAL CHARACTERISTICS OF TEST AIRPLANE 


Wing 


Total area (including flaps, slats, and 

49.92 sq ft covered by fuselage 

Aspect ratio 

Mean aerodynamic chord (wing station 98 . T in. 

measured normal to center line) 

Sweepback of 0.25-chord line 

Root airfoil section (normal to 0.25-chord line). 

Tip airfoil section (normal to 0.25-chord line) . 

. . . 287.9 sq ft 

. . . 37-1 ft 

. . . k .79 

. . . 0.51 

. . . 97.03 in. 

. . . 35°i4’ 

. . NACA 0012-64 
(modified) 
. . MCA 0011-64 
(modified) 

Horizontal Tail 


Total area (including 1.20 sq ft covered by 

Mean aerodynamic chord (horizontal- tail 

station 33-5^ in.) 

Sweephack of 0. 25-chord line 

Airfoil section (parallel to center line) . . . . 
Tail length 

. . 35.0 sq ft 

12.8 ft 
. . 4.65 

0.45 

. . 3^.7 in. 

. . 34°35’ 

. . NACA 0010-64 

18.12 ft 

Weight and corresponding center— of- gravity position 
(landing gear retracted) 


Take-off (c.g. at 23.0-percent M.A.C.) 

Landing (c.g. at 21.9-percent M.A.C.) 

Average for computation purposes (c.g. at 

22.5— percent M.A.C 

. . 14,102 lb 

. . 11,612 lb 

. . 12,800 lb 

Average moment of inertia about Y axis 

. 17,480 slug-ft 2 





Figure 1.— Photograph of test airplane showing external "boom configuration 
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Figure 2 Two - view drawing of test airplane showing research 
airspeed installation . 
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Figure 3.- Sample flight records of pressure altitude, Mach number, 
pitching velocity , and elevator angle . 
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Figure 3- Concluded . 
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Mach number, M 

Figure 4 - Variation of period of oscillation with Mach number. 
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Figure 5. Variation of damping coefficient with Mach number. 
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Figure 6 - Variation of C m ^ with Mach number. 
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Figure 7 -Variation of C m<j + C m ^ with Mach number. 
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Figure 8. - Manufacturers estimates of the variation with Mach number of lift - curve 
slope for steady flight lift coefficients . 
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Figure 9. - Number of cycles required for oscillation to damp to I/IO amplitude. 
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